Model for vortex-core tunneling spectroscopy of chiral p-wave superconductors via odd-frequency 

pairing states 
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The local density of states (LDOS) is studied theoretically in terms of the odd-frequency (odd-w) Cooper 
pairing induced around a vortex core. We find that a zero energy peak in the LDOS at the vortex center is robust 
against nonmagnetic impurities in a chiral p-wave superconductor owing to an odd-o; s-wave pair amplitude. 
We suggest how to discriminate a spin-triplet pairing symmetry and spatial chiral-domain structure by scanning 
tunneling spectroscopy via odd-oj pair amplitudes inside vortex cores. 

PACS numbers: 74.20.Rp, 74.70.Pq, 74.50.+r 



The unambiguous determination of Cooper pairing symme- 
try is of prime importance for understanding the pairing mech- 
anism of unconventional superconductivity. In the last decade, 
the discovery of the ruthenate superconductor (SC) Sr2Ru04 
has stimulated an enormous amount of studies, where a chi- 
ral superconductivity with time-reversal-symmetry breaking 
was indicated by muon-spin-rotation and polar Kerr effect ex- 
periments, and a spin-triplet state such as a p-wave one was 
suggested by Knight-shift measurements |Q],Q|. Therefore, a 
spin-triplet chiral (e.g., p x ± \p y ) state has been considered 
as a most promising pairing state [3]. However, in order to 
confirm the chiral p-wave state definitely, much more clear 
experiments based on new ideas are needed [4]. 

One way of identifying the pairing symmetry among possi- 
ble candidates is to resolve the quasiparticle states, by tun- 
neling spectroscopy, via the surface Andreev bound states 
(SABSs) H. The SABS 0] originates from a sign change 
of anisotropic pair potentials at an interface, and it is observed 
as a zero-bias conductance peak in various materials |5[]. It 
is also well known that the so-called vortex Andreev bound 
states (VABS) are formed around a vortex core [7]. These 
bound states can be interpreted as the odd-frequency (odd- 
ui) pairing states [8]. The odd-cj pairing state is character- 
ized by a pair amplitude that is an odd function of the Mat- 
subara frequency [9]. The origin of the generation of odd- 
u> pair amplitude is as follows. In inhomogeneous system, 
due to the breakdown of translational invariance, the pair po- 
tential acquires a spatial dependence which leads to coupling 
between the even- and odd-parity pairing states. The Fermi- 
Dirac statistics then dictates that the pair amplitude of op- 
posite parity should be opposite in frequency 11 1 CHI - Owing 
to an axial symmetry of a vortex, there exists a fundamental 
rule HUE] that relates the angular momentum of the odd-w 
Cooper pair and the topology (winding number) of vortex. It 
is interesting and important to confirm experimentally the ex- 
istence of odd-w pair amplitudes subjected to the topological 
symmetry rule due to the vortex. Note that the SABSs are also 
interpreted as a generation of the odd-cu pairing states iflftfol . 

In this Letter, we will show that nonmagnetic impurity scat- 



tering effect on the VABS can be used to detect the sym- 
metry rule for the induced odd-uj pair amplitudes. We will 
study the odd-cj pairing amplitude and the local density of 
states (LDOS) around a vortex for spin-singlet s-wave and 
spin-triplet chiral p-wave SCs. The effect of nonmagnetic 
impurities serves as a probe of the symmetry rule and the 
pairing symmetry in candidates for a spin-triplet chiral SC 
such as Sr2Ru04. We will also point out how to observe 
spatial chiral-domain structure by scanning tunneling spec- 
troscopy/microscopy (STS/STM). The issue of chiral domains 
has been studied experimentally in Sr2Ru04 lfl4l 0511 . How- 
ever, a more direct observation by spatially resolved probe has 
not been performed. Therefore, it is interesting to propose an 
idea to observe directly the spatial structure of chiral domains. 
In the present theory, we point out that the proposed method 
for identifying a spin-triplet chiral p-wave state can be simul- 
taneously used for detecting chiral domains. 

We consider a spin-singlet s-wave and a spin-triplet chi- 
ral p-wave pairing state with a e£-vector parallel to the z axis 
iQ 13] • In both cases, the quasiclassical Green's function 
g(iu) n ,r, fc) is represented in a 2 x 2 matrix form as [ 16] 
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It follows the Eilenberger equation 
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which is supplemented by the normalization condition g 2 + 
ff = 1. Here, u n = (2n + 1)tvT is the Matsubara fre- 



quency, vp is the Fermi velocity, and r 



,t z ) are 



Pauli matrices in the particle-hole space. We use units in 
which H = fcs = 1. r denotes the center of mass of the 
Cooper pair, and fc is the unit vector of the Fermi wave number 
(fc = fcp/lfcpl). The pair potential is expressed as A(r, fc) = 
(f x + if y )A(r,fc)/2 - (f x - if y )A*(r,fc)/2. We incorpo- 
rate the impurity scattering effect within the Born approxima- 
tion. The impurity self energy is given as iflrJl E(iw„,r) = 

ifj / )(/(iw„,r,fc)) - 



T[(r, + if„)(/(iu; n ,r,fe)) - (7 



2 



2ir z (g(iu> n , r, fc)\]/2, where the brackets (■••} denote the 
average over the Fermi surface. We define the impurity scat- 
tering rate in the normal state as T = l/2r with mean free 
path vpT. 

We assume an isotropic two-dimensional system and intro- 
duce the angle 9 as k = (cos 8, sin 9). The chiral p x ± \p y 
(— p±) wave is represented by exp(±i6>). The self-consistent 
equation for the pair potential is given by 



A(r,0) = A+(r)e +i ' e 
A±{r)=nTV 
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where V is the coupling constant (see Ref. |fl6ll for details). 
We consider a spin-singlet s-wave state (I = 0) and a spin- 
triplet chiral p-wave one (I — 1) for the pair potential. The 
Fermi-surface average is (• • • } = f Q 27r d0' ■ ■ ■ / (27r). The pair 
amplitudes represented in the Matsubara frequency are 



,r)e +iW +F^(iu; n ,r)e 
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/(ia>„,r,0')), 
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with the quantum number of the angular momentum I = 0, 
1, 2, • • • , where the even-cj and odd-cj pair amplitudes sat- 
isfy f(iu n ,r,6') = /(-iw n ,r,0')> and f(iu n ,r,9') = 
—f(—iu> n , r,8'), respectively. 

The Eilenberger equation is simplified by introducing the 
Riccati parameterization fnll . We numerically solve the re- 
sulting Riccati equations and the self-consistent equations for 
the impurity self energy and pair potential iteratively iflcill . 
Throughout the paper we set the temperature T = 0.1T C , 
where T c is the critical temperature in the absence of im- 
purities. Using the self-consistently obtained pair potential, 
we determine the self energy by the analytical continuation 
with \uj n — > E + iS 1 18]. The LDOS is then calculated as 
N(r,E) = iVp(Re g R ). Here, Np is the normal-state den- 
sity of states at the Fermi level, g K — g(\oj n — > E + iS), E 
means the quasiparticle energy, and S is an infinitesimal quan- 
tity. We select S = O.O6A0 as a typical value, where Ao is the 
bulk amplitude of the pair potential at T = and T = 0. 

An enhancement of the LDOS in the presence of odd-cj 
pairing can be understood by means of the normalization con- 
dition. Since / with E = is given by / = — /* for 
odd-cj pairing state 18J, one can show that generally N(E = 
0)/N F > 1 owing to g 2 = 1 + |/| 2 > 1. This means that 
the emergence of odd-cj pairing is a physical reason for a zero 
energy peak (ZEP) in the LDOS. 

In circular coordinates r = (r, cf>), we consider an axially- 
symmetric vortex situated at r = 0, the vorticity of which is 
perpendicular to the two-dimensional system. The axial sym- 
metry requires A — > A exp(iNa) with an integer TV for a ro- 



tational transformation 4> — * 4>+a and 6 — > 6 + a [ 1 1 
The pair potential that satisfies this constraint is 




FIG. 1 : (Color online) (a) Spatial dependence of the odd-o; pair am- 
plitudes (u) n =o), and (b) the corresponding LDOS for the s-wave 
vortex with T — 0.1 Ao (left panel) and 0.3Ao (right panel). 



for a spin-singlet s-wave vortex, and 

A(r, 6 ) = A+ ( r ) e i[e+(™-2)0] + ( r ) e i[-*+«^ , (g) 

for a spin-triplet chiral p-wave vortex. Here, m is the winding 
number of a vortex. We have assumed the p_-wave state in 
bulk for the chiral p-wave vortex, namely A+(r — > 00) = 
and A_ (7- — > 00) ^ 0. Due to the same constraint from axial 
symmetry, the pair amplitudes satisfy 



for the spin-singlet s-wave vortex, and 



(9) 



A(r,0) = A(r)e 



im(f> 



(7) 



F^(iu; n ,r) - F^Ou^r)^- 1 ^*, (10) 

for the spin-triplet chiral p-wave (p_-wave) vortex. Note that 
if a pair amplitude has a finite phase proportional to cf>, the 
amplitude inevitably becomes zero at r = where <p is unde- 
fined. 

First let us discuss the spin-singlet s-wave vortex (I = 0) 
in the presence of impurity scattering. The spatial dependen- 
cies of odd-w pair amplitudes (ui n —o) are plotted in Fig.QJa). 
There, the distance r from the vortex center is normalized by 
the coherence length £0 = ^f/Ao. In the case of winding 
number m = 1, only the p + -wave pair amplitude with odd-w 
is induced at the vortex center. It is because only the phase 
of the p+-wave pair amplitude F+~ 1 ' is zero in Eq. (O for 
m = 1. As a result, the LDOS at the vortex center has the 
ZEP as shown in Fig. [IJb), reflecting the VABS that origi- 
nates from the induced odd-w pair amplitude la [l0]. The 
magnitude of the odd-cj pair amplitude is suppressed with the 
increase of T. Accordingly, the height of the ZEP decreases 
with increasing T [Fig. [Tib)]. Actually, a collapse of the ZEP 
upon doping impurities in an s-wave SC was observed exper- 
imentally 11911 . Here, the induced odd-w pair amplitude is the 
p + -wave (i.e., non-s-wave), and therefore Anderson's theo- 
rem for non-magnetic impurities 12011 is not applicable. As a 
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result, the odd-w pair amplitude, namely the ZEP, is sensitive 
to impurity scattering. 
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FIG. 2: (Color online) Spatial dependence of the pair potentials and 
the odd-cij pair amplitudes (w n =o) for the chiral p-wave vortices with 
r = O.lAo (left panel) and 0.3Ao (right panel). The p-wave pair 
potentials are normalized by bulk value of the p_-wave one. 



Next, let us move to impurity effects on the spin-triplet chi- 
ral p-wave vortices, the main subject of this paper. There 
are two types of vortex, the vorticity of which is antiparal- 
lel (m = 1) or parallel (to = —1) to the predominant p_- 
wave pair potential 111 211 . As seen in Fig. [2] the odd-u s- and 
eL-wave pair amplitudes are induced in the vortex center for 
the antiparallel and parallel vortices, respectively. The magni- 
tudes of these odd-cj pair amplitudes are reduced with the in- 
crease of T. Note that the decrease of the pair amplitudes with 
r is weaker for the antiparallel vortex than that for the parallel 
one. As seen from Figs. [3a) and[3jb), while the height of ZEP 
strongly decreases with T for the parallel vortex, it is robust 
against T for the antiparallel vortex iflil 21]. 

The clear difference between two states can be understood 
in terms of the symmetry of the odd-w pair amplitude. Due 
to the phase factor in Eq. ( [Tol l, the odd-co pair amplitude in- 
duced at the vortex center is inevitably s-wave (I = 0) for 
the antiparallel vortex (to = 1) and <i_-wave (I — 2) for the 
parallel one (to = —1) as mentioned above. According to 
Anderson's theorem [20], an s-wave pair amplitude is robust 
against non-magnetic impurities, while a d_-wave pair am- 
plitude is sensitive to such impurities (see the inset in Fig . [5j ■ 
The VABS (corresponding to the SABS in Ref. J3) origi- 
nates from an odd-w pair amplitude 111 311 . and it is reflected 
in the ZEP [8]. Hence, the ZEP is robust against the increase 
of T only for the antiparallel vortex where the VABS orig- 
inates from the s-wave pair amplitude. On the other hand, 
for the parallel vortex the ZEP is substantially suppressed by 
r because of the d_-wave (i.e., non-s-wave) pair amplitude. 
Accordingly, measurements of the ZEP under the influence 
of impurities correspond to observations of the symmetry of 
those odd-u pair amplitudes as illustrated in Fig. [4] In order to 
visualize the striking difference between the antiparallel and 
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FIG. 3: (Color online) The LDOS at the vortex center (left) and the 
density plot of LDOS as a function of the distance r (right) for (a) the 
antiparallel and (b) the parallel vortex. The upper and lower panels 
are plots for V = O.lAo and 0.3Ao, respectively. 



parallel vortex states, we show in Fig.[5]the zero energy LDOS 
at the vortex center as a function of T. In the inset of Fig. [5] 
it is found that the odd-w s-wave pair amplitude (circle) is 
likely to follow the Abrikosov-Gorkov (AG) law [22], while 
the odd-ui d_-wave one (triangle) decays more rapidly with 
increasing T. The decrease of the odd-w s-wave component is 
involved with the fact that the odd-ui component is generated 
through the coupling to the bulk even-w p-wave pair potential 
which is sensitive to V. 

The above results have a significant implication to identi- 
fication of pairing symmetry in chiral SCs such as Sr2Ru04. 
In chiral superconductivity, the pair potential is composed of 
two degenerate components. As a result, degenerate chiral 
states, such as p x + ip y and p x — \p y , form a domain struc- 
ture in a chiral SC under field cooling condition with high- 
speed cooling rate |3- Under magnetic fields the antipar- 
allel vortex is known to be energetically favorable M23I1 . and 
therefore the antiparallel vortex state, namely one of two chi- 
ral states, dominates in a sample under field cooling condi- 
tion with slow cooling rate. In the case of high-speed cool- 
ing, a domain structure remains as a mixture of antiparallel- 
and parallel-vortex domains. If a SC is a spin-triplet chiral 
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(a) Antiparallel chiral p-wave vortex 




FIG. 4: (Color online) Schematic illustration of STS/STM measure- 
ments to detect the odd-o; pair amplitude for (a) the antiparallel and 
(b) the parallel vortex in chiral p-wave SCs. The arrows represent the 
phase rotation. 
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FIG. 5: (Color online) The ZEP in the LDOS at the vortex center as 
a function of V. The lines are guides for the eye. Inset: the odd- 
ui pair amplitudes (w n =o = ~kT) at the vortex center as functions 
of r. The s-wave (circle) and d--wave (triangle) pair amplitudes 
are plotted for the antiparallel and parallel vortices, respectively (see 
also Fig. [2}. Each pair amplitude is normalized by its value at F — 0. 
The solid curve represents the Abrikosov-Gorkov law l22ll . 
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p-wave one, such a difference in the chiral state between slow 
and high-speed coolings is observable via the ZEP at vortex 
cores shown in Fig. [4] In addition, through a distribution of 
those two different vortex states, the existence of a chiral do- 
main structure can be observed by STM at zero-bias under 
high-speed cooling with a spatial resolution of the order of 
inter-vortex distance. On the other hand, for any spin-singlet 
pairing state such as a chiral <i-wave one [e.g., p z (p x ± ipy)], 
the induced odd-w pair amplitudes are inevitably spin-singlet 
odd-parity ones [8] and the odd-w s-wave (= even-parity) pair 
amplitude is never induced. In Sr2Ru04, the ZEP at a vortex 
core has been observed B24I1 . and samples with different impu- 
rity scattering rate could be prepared B22I1 . Therefore, an ex- 
perimental setup proposed here may provide strong evidence 
for spin-triplet chiral p-wave superconductivity. 

In conclusion, we have studied the odd-cj pair amplitudes 
and the LDOS around a single vortex in s-wave and chiral 
p-wave SCs. For the antiparallel and parallel chiral p-wave 
vortices, we have found that the odd-w s-wave and <i_-wave 
pair amplitudes are inevitably induced at the vortex center, 
respectively. The robustness of the ZEP at the vortex core 
against non-magnetic impurities originates from this odd-uj 
s-wave pair amplitude. Those odd-ui pair amplitudes can be 
observed by STS/STM, serving as a probe of a spin-triplet 
pairing and of a spatial distribution of the chiral domains. 
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